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Abstract 

In this paper we discuss a model which can be considered as a generalization of the 
well-known scalar-tensor Brans-Dicke theory. This model possesses an interesting feature: 
due to Born-Infeld type non-linearity of the scalar field the properties of the interaction 
between two test bodies depend significantly on their masses. It is shown that the model 
| can be interesting in view of the Pioneer 10, 11 spacecraft anomaly. 

o 

o ; 1 Introduction and setup 

One of the most known scalar-tensor theories of gravity is the Brans-Dicke theory [1, 2j. It 
^2 • describes the scalar field non-minimally coupled to gravity with the action 
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(3JT)! where uj is the Brans-Dicke parameter and L ma u er is the Lagrangian of matter. In the limit 
*~ ~ '■ uj —>■ oo the theory goes to the standard General Relativity. This theory is very well examined, 
the present days gravitational experiments set stringent limits on possible values of the Brans- 
Dicke parameter uj [3l H] . 

■^j- \ In this paper we consider a generalization of this theory based on the use of the Born-Infeld 
scalar field. This field itself was widely discussed in the literature, see, for example, [5] [TT] 
and references therein. We will show that such a highly non-linear covariant theory possesses 
very interesting features and, in principle, it can account for the anomalous acceleration of 
Pioneer 10 and Pioneer 11 spacecraft [121 [32], leaving the planets of the Solar System devoid 
of such extra constant acceleration which is excluded by observations [13] ■ 

To this end let us consider the following four- dimensional action describing Born-Infeld 
scalar field non-minimally coupled to gravity 



S = d^xyf^g 



, uj g^d^d v y 

ipK + J \ / 1 J + Lmatter 



f V 

with / > 0. The action of the scalar field has a non-standard form, but in the limit / — > oo 



(2) 



„ uj g^d^ipdutp uj g^dpydvtp 

hll ~l~ yf ~2—^' 

and (T5]) transforms into the well-know action of the Brans-Dicke theory ([T]) with the Brans-Dicke 
parameter uj = ^. The theory with action (T5]) can be considered as a generalized Brans-Dicke 
theory. The extra term — / in (j2J) is added to preserve the Minkowski background metric in the 
vacuum state. 
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The vacuum expectation value of the field ip is supposed to be (p vac = M Pl . Let us represent 
the scalar and the gravitational fields as 



V = M Pl + ^0, (3) 



g^v = + -tjt V; ( 4 ) 

Mpi 

where mc = 0, g^f = rj^, r]^ u = diag(—l, 1, 1, 1) is the flat Minkowski metric, and expand 
action (j2J) into series with respect to <fi and h^. The "gravitational" part of the action 



d A x\/—gtpR 



can be represented as 
d A x 



M Pl L^[h, v ] + -L</>L«[V] + ( 1 + T7^0 ) £ (2) [V] + - 
^ V m pi v u 



(5) 

where l/ 1 ^/^] is linear in h^ u , [h^] is quadratic in h^ v and so on. The Born-Infeld part of 
the action 

, 4 , — f L w g^d^dyip 



d xy/=gf\ 1 



f V 



can be rewritten as 



<;**/,/ 1 - j (>r - jj^w + ...) b^4 (i - ), (6) 



where we have omitted the expansion of sj—g. Here and below indices are raised by rf" . 

Now let us discuss formulas (j^J) and (jSJ). First, the term L^[h^ v ] is simply L^/t^] = 
d tJl d v h IM , — d^d^h (where /i = /i^), which is a total derivative. Thus, the term M P iL^[h^ u ] 
vanishes from the action. Second, since we suppose to work in the Newtonian approximation, we 
can drop the term j^-h^" (as well as higher corrections in h^ u ) in comparison with rf v . For these 

reasons we can also drop the terms Z/ n ) [/t^] for n > 2. As for the term A/ ^^ 0, it is not evident 
that it is much smaller than unity. Nevertheless, let us suppose that j^— j^ (j> "C 1 and drop 
the term jj^-^4> and the subsequent terms in (jSJ), as well as the cubic term M ^^ <pL ( ^ [h^ u ] . 
Below we will show that condition ., 1 ^ <b <C 1 indeed holds. Note that we are not able to 

drop the quadratic term ^=Z/ 1 < ) [h^ u ] because it ensures the interaction of Born-Infeld scalar 
field with matter. 
Thus we get 



Seff = j d 4 x Uj=LV[h^\ + L^ih^] + f^j\--^d^d v <P+ (7) 



1 -hTt, 



2Mf " ' ' 
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where £„„ is the energy-momentum tensor of matter and 



L^[h flu ]=L FP [h fiU ] = (8) 



™ [d P h^d p W v - d p hd p h + 2d^ v d u h - 2d^ v d p h pv ] 
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is the standard Fierz-Pauli Lagrangian. It is convenient to diagonalize action (JTj) with the help 
of the standard redefinition 

h/u, = V j=V^- (9) 

After some algebra we get 



L eff = L FP [b^ v ] + f x /i - j-q^d^d^ - ^-^d^d^ + (io) 

1 -\rt 1 



2Mpi 2Mpi^ ' 

where t = rj^ u t^ u . The extra term — -^rj^d^dycf) in (J2J) has appeared in the action after 
diagonalization. 

The Born-Infeld part of action (fTUj) has the standard form of the Dirac-Born-Infeld (DBI) 
scalar field action, though the standard DBI Lagrangian has a different origin. It is necessary 
to note that we take / > (like in [8]), contrary to the case / < 0, which is often discussed in 
the literature (see, for example, [5]-[7], [TO]). 

It should be also noted that we neglected the term — / (see action (jSJ)) while obtaining 
Lagrangian fllOl) . We will discuss this issue in the next section. 

Lagrangian (1101) allows one to examine the stability of the model at least above the Minkowski 
background. To this end we consider t^ v = and suppose that rf v d^d v <^) <C /. Expanding 
the square root in ffTOj) up to the linear term we get a quadratic Lagrangian 

L 2 = L FP [b, u } - — — rTd^t. (11) 

For ^> 1 (this case will be discussed below) the kinetic term of the scalar field has the proper 
sign, which leads to the absence of ghosts in the theory. Indeed, for rj^d^d^ <S / Lagrangian 
(II Op describes the standard Brans-Dicke theory in the Newtonian approximation, which is 
known to be stable. Higher corrections to (1111) lead to an infinite tower of self-interaction 
and interaction terms. Nevertheless, the solution which will be discussed below corresponds 
to a deep non-linear regime of the Born-Infeld part of the model where perturbation theory 
does not work. It has appeared to be very difficult (maybe even impossible) to examine the 
perturbations around this solution analytically, one should make numerical analysis. Thus, the 
question about the stability of the solution presented below has no definite answer yet. 

2 Equations of motion and extra anomalous force 

The equations of motion following from Lagrangian (fTU|) take the form 

U K» = (v - \v»ut) , (12) 




^ .V<W = ' «, (13) 



1 - ^d p (j)d a (f) w / 2M piV^ 



f 

where □ = rj^d^dy. We used de Donder gauge condition d^b^ u — \d v b = while obtaining 
(ED. 
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We will be interested in the additional interaction between two bodies caused by the scalar 
field cf). As can be seen from initial action (j2j), ordinary matter interacts only with the metric, 
in this sense the weak equivalence principle is fulfilled. Thus the gravitational force acting on 
a test body can be easily obtained by considering geodesic motion, and in the Newtonian limit 
we get the well-known formula 

^2^ VA »- (14) 

If one considers a non-point-like source, this formula transforms into a formula describing the 
force acting on the center of mass of the test body 

mR = — — I dVp{x)Vh 0Q , (15) 



2M Pl Jy 

where p(x) is the density of the body of volume V such that 

p(x)dV = m 

v 

and R is the vector pointing to the center of mass of the body. Using (JHJ) we get the standard 
Newtonian force (the contribution of b^), and an anomalous force 

dV>(£)V0 (16) 



2M Pl ^J Jy 



(the contribution of <fi). We will be interested in this anomalous extra force. 

Now let us turn to Eq. (|T3i) and consider the static case of spherically symmetric bodies. 
We suppose that the energy-momentum tensors of the test bodies have the form 



too = pl,2(x), <i/ = 0- 



Due to the spherical symmetry 



p 2 {x) = p(r), r < r*, (17) 
p 2 (x) = 0, r > r„, (18) 

where r* is the radius of the second body. 

Let us denote rj^djtp = V0 = <p. Then Eq. (1131) for the case of two test bodies takes the 
form 

1 " 4tt (- j_ Pl (x) - X p 2 {x) ) . (19) 

Now we are ready to examine the force acting on the test body in such a system. The 
coordinate system that will be used for calculations is presented in Fig. [TJ The force will be 
calculated for the second body (the right body in Fig. [Q . 

The solution to Eq. (fT9|) inside the second body has the form 

^ 3 r 1 (Ml m(r)f . 

+ -^ = ~ o_ /=I7 hr + -T- ■ (2°) 
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m 




where I = R + r, r = M is the mass of the first body, 

pr 

m(r) =An f 2 p(f)df 



for r < r* and 7n(r)| r > n = m (see Fig. [T]). 

It is convenient to represent the parameter / as 



/ 



M 2 



(21) 



where Tm is a parameter depending on M. 

We suppose that a; ^> 1. Since y = j= > 1 and 3/u <C 1, we can look for a solution to 



Eq. (I20I) using the perturbative approach. After some algebra one can get 



V7lf + 



Ml I m(r)f 



M 2 , I Ml , m(r)f 

7£ + [ + ,.3 



i-i 



M 2 



M 2 , [ Ml"-)- m(r)f 



/ 



:22i 



For the objects which will be discussed in the next section of the paper the correction in (1221) 
~ 3/u appears to be at least ~ 10~ 3 or even smaller. Thus we can drop this correction and use 



Ml m(r)r 
1 — 

/3 T r 3 



v7 



(23) 



M 2 

^17 



Ml _|_ m(r)f 



For tm 3> -R (1231) transforms into 



Ml _i_ rra(r)r 
23 "T" r 3 



M[* I m(r)f 
;3 t" r 3 



(24) 



We can estimate the maximal value of the field (f> itself. The approximate size of the non- 
linearity zone is Tm, and (ft can be estimated as 



M 



87T^/LjMpir M 



and thus 

M 

n ^ J . 

A more accurate analysis based on the use of the solutions inside and outside the non-linearity 
zone (the latter behaves as ~ 1/L, where L 3> tm is a characteristic distance from both bodies), 
provides an analogous estimate (up to the factor of the order of unity). For the parameters, 
which will be considered in the next section, ^f Mpi "C 1 and the corresponding terms in (151) 
indeed can be omitted. 

Now let us estimate the effects that could be produced by the omitted term — / of action 

(J2J) (see previous section). In the non-linearity zone y^l — j((f(f)) <C 1 and the term — / is not 
compensated. It indicates that the background metric in the non-linearity zone is not the flat 
Minkowski metric, but a de Sitter-like background metric leading to a local expansion. For the 
observer, say, on the first body it looks like a repulsive force acting on the second body. This 
force has the form 

f - 

F T ep ~ m—p^-R, (25) 
M Pl 

which can be easily obtained by considering de Sitter metric in the static form [2]. Using (T2~TT) 

we get 

/ M \ ( R\ 1 M 



- 1 rep 



( — ) (*.) 

\^M 2 pi r M ) \r M J 



16ttM pi ur 2 M 



m. (26) 



For the values of the parameters that will be used below this force appears to be much smaller 
than the forces caused by the fields b^ v and obeying (TT2"|) and (IT31 respectively. Thus, for our 
purposes we can use the Minkowski background metric instead of a de Sitter-like background 
metric. We would like to note that analogous estimates can be obtained if we retain the term 
—/in the action and get slightly modified equations for the fields b^ u and </>. 



3 Specific examples 

Now we turn to the effects which can be produced by the DBI scalar field in our Solar System. 
Let us suppose that 

M = M , r M « 100 AU, u w 700, 

which means that / « 2 • 10" 44 GeV 4 (our "reduced" Planck mass M Pl ta L2 - 1 °^g eV « 1.7 ■ 
10 18 GeV). 

In what follows we will consider two cases: 

1. A light body with the mass m, for which ^ ~ ^±r- <C jp, r* <C R <C tm (for example, 
a spacecraft like Pioneer 10, 11 with m ~ 300%, r* ~ lm). In this case in the leading 
order 



f Ml I m(r)f \ ' 
\ «3 "+" r-3 y 



MZ m(r)f \ v7 /! f ^ rn(r)rR 2 

7^ + r 3 j "7(Mp = \i? + Mr 3 
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It is worth mentioning that there is no such static solution for the case / < 0. Indeed, if 
/ < then r\j < and we get negative values under the square root (see Eq. fl23J)), which 
is the consequence of the existence of a horizon at a finite distance (see, for example, [9], 
where solutions with horizons in DBI scalar field theory are discussed). That is why we 
chose the case / > 0. 



Substituting the latter formula into (1161) and integrating over the volume of the body 
leads to (we use the fact that J rdfl = 0, where fl is the solid angle) 



1 



M 



IQtiM 2 pi ujr 



m. 



(2? 



It is necessary to note that in principle 



. dVL ^ 0. 



(29) 



Ml , 
i 3 r 



m(r)r 
3 



But we can neglect possible corrections because anyway 

Ml 



I 3 



> 



m r r 



see 021]). 
Formula 



is written in the system of units H — c — 1. The replacement 



G 



allows one to pass to the SI units, which results in the acceleration towards the Sun in 
the SI units 



GM 



8.7 • 10- 10 m/s 2 . 



uor 



(30) 



M 



We note that this acceleration does not depend on the distance from the Sun, which is 
exactly the situation with the Pioneer 10 and Pioneer 11 spacecraft [T2j [T3"]. 

As for the bodies on the surface of the Earth, we can carry out analogous calculations 
taking M = M© (in this case Tm also changes). Our ideal test bodies with density 
profile (fTTl) . ffTBl on the surface of the Earth also possess an additional acceleration a anom 
towards the center of the planet. It is evident that this acceleration can be neglected in 
comparison with g w 9.8m/ s 2 for Earth-based gravitational experiments. 



2. Heavy bodies with the mass m (planets), 



m(r) 



<C R tm- In this case we 



should carry out calculations more precisely because 



Ml 



m(r)f 



and possible corrections due to (j29j) can be quite large. Correspondingly, we should take 

Mr(fR) 



\ 



Ml ///( /•)/•' \ m [ r) 
~P~ + r 3 



1 + 2 



Mr(fR) m(r) 
m(r)R 3 r 2 



1 + 



m{r)R? 



(31) 
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and thus 



m (r r 



v7 



Finally we obtain 



Ml I m(r)f 
,3 "T _3 



MZ m(r)r 
1 — 



/I- 



Mr 5 



m(r)R 2 r \Rr 



(32) 



>/7 



+ 



Mr 5 



r m(r)R 2 



Ml , m(r)f 

is -r r 3 



R f I Rf 
R r \ Rr 



(33) 



Substituting the latter formula into (fIBl) and integrating over the volume of the body 
leads to 



1 



8vrM 



167rMp Z 3ur M m 



mir 



MmR 
R 3 ' 



in the SI units 



where 



-G 



MmR 



eff 



R 3 



G 



eff 



G 



8vrM 
2>uir 2 M m 



'• p(r) 



r dr. 



(34) 

(35) 
(36) 



One can see that the extra force acting on a heavy body ~ 1/R 2 . Such a behavior is 
inherent to the ordinary Brans-Dicke theory and we can replace the original potential 
<p in (fTB"]) by an effective potential ~ 1/R. The effective Brans-Dicke parameter can be 
easily extracted from (136]) : 



1 



8ttM 



2ujbd + 3 Zuir 2 M m J m(r 



r dr, 



^BD 



2>ujr 2 M m 



16ttM f n r * 4\r 4 dr' 

JO m(r') 



(37) 



(38) 



A significant difference from the original Brans-Dicke theory is that ujbd depends on the 
mass m, i.e. it is different for different planets. 

To estimate ujbd for different planets we suppose that p(r) = t=3? = const. In this case 



'Lwri 



umr 



^BD 



M 



a ff 



2Mrl 



2a n 



(39) 



where a// is the free fall acceleration on the surface of a body (a planet). For example, 

l9 



^Mercury « 2 .1 • 10 9 



Jupiter 
BD 



w 1.4- 10 10 . 



Such large values of the Brans-Dicke parameter do not contradict the experimental bounds 
ujbd > 3500 obtained in the Solar System gravitational tests [31 H] (we would like to note 
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that the bounds on the ~ R/R 3 extra force differ from the bounds on the ~ R/R extra 
force) . 

For the other limiting case p(r) = 2 ™2 r 

co BD « (40) 

One should note that in the case of a heavy body with p(r) = const there exists a region 
r < f such that ~ For this region one should carry out the calculations described 
in item 1 (the case of a light body). But even for Mercury (if we suppose p(r) = const) 
y- pd 1.8 ■ 10~ 2 , and the constant extra acceleration appears to be 

f 3 _ 6 

O'extra ~ "J O'anom ~ 6 • 10 Q>anomi 

which does not contradict the existing experimental restriction on a possible extra con- 
stant acceleration of the planet [H] (which should be much smaller than that of the 
Pioneers 10, 11 spacecraft). It is easy to check that for the other planets of the Solar 
System a exira also does not exceed the experimentally allowed limits [13] . If one takes 
p ~ - this region is absent and thus a extr a = 0. Of course, our density profiles for the 
planets are an idealization, but more realistic profiles should lead to the values of a extra 
which lie somewhere between the limiting values obtained above. The same is valid for 
the values of the effective Brans-Dicke parameter lobd- 

The physical difference between the two cases can be easily explained. In the first case vectors 
</> are directed approximately parallel to the vector R at any point of the second body, whereas 
in the second case these vectors are approximately parallel to the radius-vectors r, which leads 
to the result discussed above. 

Thus we have shown that in principle it is possible to construct a covariant theory which 
"distinguishes" light and heavy test bodies with respect to an external gravitational field of 
a source. We would also like to note that the value of the parameter / which is necessary 
to reproduce the anomalous acceleration of the Pioneer spacecraft was chosen to be / ~ 2 • 
10~ 44 Ge^ 4 . This value is quite close to the vacuum energy density ~ 10~ 47 GeV 4 , responsible 
for the accelerating expansion of the Universe. In this connection it is very interesting to 
examine possible cosmological manifestations of the model described by action (T5]). This issue 
calls for a more detailed and thorough investigation. 

Of course we can not argue that the Pioneer anomaly is indeed caused by the existence 
of such a DBI scalar field. Moreover, recently it was shown that a part of the anomalous 
acceleration can be explained by the thermal recoil force effect [15]. Nevertheless the model 
presented in this paper possesses quite interesting features, does not contradict experimental 
data at least in the Newtonian limit and seems to be worth an additional examination. 
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